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ABSTRACT 

We present observations of ro- vibrational OH and CO emission from the Herbig Be star HD 100546. 
The emission from both moleeules arises from the inner region of the disk extending from approx- 
imately 13 AU from the central star. The velocity profiles of the OH lines are narrower than the 
velocity profile of the [O I] A6300A line indicating that the OH in the disk is not cospatial with the 
O L This suggests that the inner optically thin region of the disk is largely devoid of molecular gas. 
Unlike the ro-vibrational CO emission lines, the OH lines are highly asymmetric. We show that the 
average CO and average OH line profiles can be fit with a model of a disk comprised of an eccentric 
inner wall and a circular outer disk. In this model, the vast majority of the OH flux (75%) originates 
from the inner wall, while the vast majority of the CO flux (65%) originates on the surface of the 
disk at radii greater than 13 AU. Eccentric inner disks are predicted by hydrodynamic simulations of 
circumstellar disks containing an embedded giant planet. We discuss the implications of such a disk 
geometry in light of models of planet disk tidal interactions and propose alternate explanations for 
the origin of the asymmetry. 

Subject headings: stars: circumstellar matter, individual (HD 100546), protoplanetary disks 



L INTRODUCTION 



An important tool for tracing the evolution of circum- 
stellar disks is the measurement of t he spectral energy 
distribution (SEP) of young stars (e.g. Strom et al.lll989l : 

19981). The differ- 



iHillenbrand et al.lll992l : iMalfait et al.l 



ences in the SEDs of young stars have led some to propose 
that they reflect an evolutionary sequence from optically 
thick disks to transitional disks (i.e. disks with optically 
thick outer disks and optically thin inner disks) to opti- 
cally thin disks. A popular interpretation of transitional 
objects is that they reflect ongoing planet formation, but 
there are other processes that can also give rise to disks 
with transitional SEDs. 

The proposed scenarios for the origin of transitional 
disks can be dis tinguished by the r adial distribution of 
gas m the disk (jNaiita et al.ll2007[) . The grain-growth 
and planetesimal formation scen arios predict an o pti- 
cally thin, gas-rich inner disk fe.g. lStrom eTaI|[T989t) . A 
system that is forming a Jovian mass planet will show 
a radial gap in the gaseo us inner disk in the v icinity 
of the planet's orbit (e.g. ISkrutskie et al. Ill990f ). Sys- 
tems that have formed massive companions (>5 Mj; e.g. 



iLubow et all 119991) or ar e being photoevaporated away 
(e.g. iClarke et al. |[200l will be depleted in gas in the 
optically thin region of the disk. Thus, we can deter- 
mine which scenario likely applies to any given transi- 
tion object by probing the distribution of gas in the in- 
ner disk. As the planet-forming region of disks is gen- 
erally spatially unresolved, high-resolution spectroscopy 
becomes a surrogate by spectrally resolving the velocity 
of gas. Assuming Keplerian rotation, spatial information 
can be extracted if the stellar mass and disk inclination 
are known. 

Work to date that takes this approach presents con- 
flicting results for well-known intermediate-mass tran- 
sition objects such as HD 100546. HD 100546 is a 
nearby (d=103t6 Pc; Ivan den Ancker et al.lll997h Her- 
big Be star with an ostensibly massive disk (Mdisk = 



(jBouwman et al. 


1 i\ — 

20031) and imagi 


Space Telescope ( 


Gradv et al.ll2005D 



indicate the presence 
of an inner hole in the dust distribution at 10-13 AU. 
While ro-vibrational CO observations of HD 100546 
show that there is no CO within 13 AU of the star 



(|Brittain et al.ll2009l: Fvan der Plas et al.ll2009D . evidence 
of ongoing accretion indicates som e gas persists in the 
inner hole (|Guimaraes et al.ll2006[) . 

In addition, the spectral profile of the [O I] A6300A 
line indicates that some gas extends inw ard to at 
least 0.8 AU (lAcke fc van den Ancke^ I2OO60. Follow- 



ing the work by iStorzer fc HoUenbachl (I1998I. 2000), 



lAcke fc van den Anckerl 1 2006f ) suggest that the O I 



emission may arise from the photodissociation of OH 
If this is the case, then the inner 13 AU of the disk is cu- 
riously rich in OH while depleted in CO, cont rary to pre- 
dictio ns for the chemistry of inner disks (e.g. iThi fc BikI 
[20051) . 

Rotational OH emission has been detected from the 
disk surrounding HD 100546, however, the spatial and 
spectral resolution are too low to resolve the line and 
thus pinpoint the location of the gas within the disk 
(jSturm et al.l I2010D . Higher excitation ro- vibrational 
OH emission has been reported from other young stars: 
from the classical T Tauri stars AS 205A and DR. Tau 
(ISalvk et al.| [2008l) as well as SVS 13jNaiita et al.l[200l 
and V1331 Cyg (jNajita et al.l 120071 : |Doppman n et af 



2011 ) , and several Herbig Ae/Be stars (iMandcU et ah. 
20Cia iFedele et al.|[20Tl[) . iMandell et al.l (|2008[ ) report 



P2.5 through P9.5 ro- vibrational lines of OH for two ob- 
jects: AB Aur and MWC 758. 

In this paper we report the observation of ro- 
vibrational OH and CO emission from the transitional 
disk surrounding HD 100546 (Sec. [3]). We compare 
the line profile of the OH emission to that of the ro- 
vibrational CO emission and the [O I] A6300A emission 
line in order to explore the properties in the region of 
the disk from which these emission features arise (Sec. 
0]). A detailed comparison of multi-epoch CO observa- 
tions of this source is presented in a companion paper 
(|Brittain et al.ll2012D (hereafter Paper II). 

2. OBSERVATIONS 

High-resolution (R=50,000), near- infrared spectra of 
HD 100546 were acquired December 22, 2010 and De- 
cember 23, 2010 using PHOENIX at the Gemini South 
telescope with th e four pixel (0".34) slit (jHinkle et al.l 
[20031 [20001 [l998h . A summary of observations is pre- 
sented in Table [TJ The position angle of the slit was 90° 
East of North for all observations. The spectra are cen- 
tered at 3145 cm-i, 2844 cm~\ and 2032 cm"! (Figs. 
(Hiaill). The PSF of the continuum is 0."7. The spectral 
settings were selected to minimize telluric absorption of 
the CO and OH emission lines. The M— band observa- 
tion includes the v=l-0 P25, P26, and P27 lines of CO 
as well as numero us other hot ba nd transitions (A v=l, 
v' > 2). Sec Britta in et al.l (120091 ) for complete line iden- 
tifications. The i— band observations cover the P9.52±, 
P10.5i±, P15.52±, and P16.5i± v=l-0 lines of OH. The 
notation is described in the appendix. 

Because of the thermal background in the L— and 
M— bands, the data were observed in an ABBA se- 
quence where the telescope was nodded 5" between the 
A and B positions. The observations were combined as 
{A — B — B + A)/2 in order to cancel the sky emission to 
first order. Each frame was flat fielded and scrubbed for 
hot pixels and cosmic ray hits. The sky lines were tilted 
by about a pixel along the full length of the slit. The data 
were rectified in the spatial direction by fitting a polyno- 



mial to the centroid of the point spread function (PSF) 
of the continuum. The data were rectified in the spec- 
tral direction by fitting a sky emission model to each row 
and interpolating the wavelength solution of each row to 
the row in the middle of the detector. To calibrate the 
spectrum, we fit a telluric model to the spectrum. The 
model was generated by t he Spectral Synthesis Program 
()Kunde fc Maguird [TQTI . which accesses the 2 000HI- 
TRAN molecular database (jRothman et al.ll2003[ ). The 
fit to strong lines is generally accurate to within a frac- 
tion of a pixel (~1 km s~^). However, near saturated 
lines (particularly water lines), there is a degeneracy be- 
tween the width of the line and the dispersion; the latter 
varies across the array. The level of uncertainty in the 
calibration near such saturated water lines is typically at 
the level of a pixel (^^1.5 km s~^ for Phoenix). Of the 
four OH lines we observed, one line, P10.5i+, was on the 
wing of a broad, saturated telluric water line. We shifted 
this line by 1 pixel blueward relative to the other three 
lines such that the standard deviation of the residuals 
between this line and the other three OH lines corre- 
sponded to the standard deviation of the noise along the 
continuum. 

We extracted the one dimensional spectrum of the star 
weighted by the stellar PSF. The extraction window was 
13 rows (I'M). The spectrum of HD 100546 was ratioed 
with the spectrum of the standard star. Regions of the 
spectrum where the transmittance falls below 50% are 
omitted (Figs. [Il[2l|3l). The equivalent widths of the OH 
lines arc reported in Table 2, and a summary of the CO 
line positions is reported in Table 3. 

3. RESULTS 

3.1. OH Temperature 

We detect the v:=l-0 P9.52±, P10.5i± OH lines (Fig. 
O and not the v=l-0 15.52±, P16.5i± lines (Fig. ^. We 
use the higher rotational lines to place an upper limit 
on the rotational temperature of the gas. If the gas is 
optically thin, the flux of an emission line is related to the 
column density of molecules in the upper energy state. 



F 



hcujjAijNj 
47rd2 



(1) 



where Ni is the number of OH molecules in the upper 
energy state, hci^ij is the energy of the emitted photon, 
and Aij is the Einstein A coeiEcient. For gas in local 
thermodynamic equilibrium, the relative population of 
the energy levels is given by the Boltzmann distribution. 



N, = Ng,e 



IQ 



(2) 



where gi is the statistical weight of the upper energy 
state, Ei is the energy of the upper state, N is the to- 
tal number of molecules, and Q is the partition function. 
Thus ln{Fij/hcf~ijAijgi) = -E,/kT + ln{N/An(fQ) so 
that the rotational temperature is given by the nega- 
tive reciprocal of the slope of a linear least squares fit of 
ln{Fij / hcUijAijg.i) and Ei/kT. Since we only detect the 
P9.52±, P10.5i± lines we can only place an upper limit 
on the rotational temperature of the gas. We fit the de- 
tected lines and the la upper limit of the non-detected 
hues and find that T^t <1400K (Fig. ^. This upper 
limit is consistent with the range of CO rotational tem- 
peratures at the inner edge of the disk (T = 1400 ±400 K; 



3 



iBrittain et al.|[2"009[ ). As the critical densities for rota- 
tional excitation of the levels we observe are high, it 
is possible that the OH is not rotationally thermalized. 
However, it seems likely that a disk environment that 
is dense enough to thermally populate the v = 1 vibra- 
tional state of CO (jBrittain ct al.ii2009h would be dense 
enough to populate the rotational states of OH. 

3.2. Comparison of HWZI of CO, OH, and 01 

We previously reported the January 2006 detec- 
tion of ro-vibration al CO emission from HD 100546 
(|Brittain et al.ll2009D . In our earlier report, we showed 
that the gas is excited by a combination of UV fluo- 
rescence and coUisional excitation and that the emitting 
region extends from 13 AU to at least 100 AU. In data ac- 
quired five years later, which are presented here, we find 
that the equivalent width of the fundamental CO emis- 
sion lines increased by 50% whereas the average equiva- 
lent width of the hot band lines increased by 30% (Paper 
II) . We also find that while the line profile of the hot band 
CO emission lines remained constant, the line profile of 
the v=l-0 emission lines varied. In the present epoch, 
the v=l-0 CO P26 line profile is symmetric, whereas the 
majority of the luminosity in the hot band CO lines is 
on the blue side of the lines. 

A detailed study of multi-epoch CO observations is 
presented in Paper II. In that paper we show that over 
three epochs spanning 2003-2010 the spectral line pro- 
file and equivalent widths of the v=l-0 CO lines vary; 
the equivalent widths of the hot band lines are less vari- 
able and their spectral line profile remains constant. Be- 
cause variations in the M-band continuum can account 
for the modest variations in the equivalent width of the 
hot band lines, we assume that the hot band lines are 
approximately constant in flux, and that the v=l-0 lines 
vary more dramatically. In this paper we restrict our 
comparison of the OH and CO line profiles to data taken 
contemporaneously, and we focus on the hot band CO 
lines because they are less variable. 

Because the individual OH lines are of low signal to 
noise (a result of their low equivalent width) we con- 
structed an average profile of the unblended OH lines in 
order to facilitate the comparison of the OH, CO, and 
[O I] A6300A fine profiles (Fig. [5] and Fig. ^. 

The fuU width at zero intensity (FWZI) of the OH 
emission is 24 km s~^ (Fig. [7|) and is consistent with 
the CO emission (Fig. [8]). The [O I] A6300A emission 
line (Fig. [T]) has a FWZI of 90 km s~|. Wh ile the O I 
extends to about 0.8 AU ()Acke et al.ll2005[) . the inner 
extent of the OH is consist ent with th at of the CO which 
is truncated at 13 ±6 AU (IBrittain e t al. 2009). 

The profiles of the OH lines differ from the profiles of 
the hot band CO lines. Inspection of the average CO 
line profile reveals a modest asymmetry where the emis- 
sion peaks on red side of the line, though the blue side 
contains the majority of the flux. This is s imilar to what 
was observed in 2006 (jBrittain et al.ll2009f ). The OH line 
profile is also resolved and asymmetric (Figs. [SI [71 [H])- 
The flux ratio of the blue side of the OH line to the red 
side of the OH line is 

We determine the Doppler shift of these asymmetric 
lines by fitting our model to the data. The barycentric 
Doppler shift is -1-16 km s^^. This is similar to the radial 



velocity we in fer from submillim eter observations of CO 
(~14 km s~^: [Panic et al.ll2010f ). It is also similar to the 
value inferred by observation of the [O I] A6300A line 
(-18 km s-^: IAcke et all [20051) . 

3.3. Origins of the Asymmetry 

We consider four possibilities for the origin of the asym- 
metry seen in the average OH line. 

Firstly, a wind may gener ate the blue-shifted asym - 
metry seen in the OH line. IPascucci fc Sterzi3 ([2009D 
show that the [Nc II] 12.81 /im emission in multiple T 
Tauri stars has a blue-shifted centroid of approximately a 
few km s~^ and a slight blue excess that they attribute to 
a disk wind. One difference between the [Ne II] emission 
and the OH line profile studied here is that the [Ne II] 
line profiles are marginally asymmetric, whereas our OH 
line profile is dramatically so. 

Molecular line profiles indicative of disk winds have 
also been reported, although the line profiles are also 
much more symmetric than our OH line profile. In their 
study of CO emission from T Tauri stars. Bast ct al. 
(2011) reported symmetric CO line profiles with a nar- 
row peak and broad wings for a subset of the sources they 
studied. The spatial constraints on the emission led the 
authors to conclude that the emission likely arises from 
a combination of emission from a circumstellar disk and 
a slow-moving disk wind. Detailed modeling by Pontop- 
pidan et al. (2011) shows that such a combination can 
produce a reasonable fit to the observed line profiles and 
their spectroastrometric properties. Additionally, NIR 
emission from organic molecules has been detected for 
several T Tauri stars, with line profiles similar to that 
found for the CO emission, suggesting a similar origin 
(Mandell et al. 2012). For systems with large disk in- 
clinations like that in HD 100546, self- absorption is ex- 
pected if the wind component is strong, which is not seen 
in HD 100546 (i=45°). Ro-vibrational CO emission lines 
forming in outflowing material of V1647 Ori have also 
been observed, but these have a P Cygni profile an d a 
symmetric emission component ([Brittain et al.l [20071 ). a 
different line profile from that of the OH in HD 100546. 
While we cannot rule out that the ro-vibrational OH lines 
we observe form in a disk wind that gives rise to a blue- 
shifted excess, our profile is distinct from the other pro- 
files in the literature attributed to a disk wind. 

Secondly, transonic turbulence in the disk (e.g., driven 
by the magnetorotational instability) produces spatial in- 
homogencities in the temperature and density distribu- 
tion of the disk. If they are large, such inhomogeneities 
might produce an OH li ne profile as asym metric as ob- 
served. Simulations by iFromang fc Nelso n (2006) find 
la fluctuations from the mean density of —10%. A much 
larger fluctuation is likely necessary to produce the level 
of asymmetry that we observe. An additional potential 
difficulty with this picture is that we observe a much 
more dramatic asymmetry in the OH line profile than 
the CO profile. Could both profiles arise from the same 
turbulent disk? If the critical densities of CO and OH 
for ro-vibrational excitation differ significantly, their sen- 
sitivity to inhomogeneities might be more pronounced 
for one molecule than the other. This picture could be 
tested with future observations. If the asymmetry orig- 
inates from turbulence in the disk, the OH line profile 



4 



would be expected to show stochastic hne profile varia- 
tions, and we would not continue to see the profile that 
we now observe. 

Thirdly, a localized hot-spot or a dense lump of ma- 
terial could induce non-axisymmetric emission. Such a 
hot-spot in the disk could be caused by a planet/disk in- 
teraction if the planet warms the inner rim in its vicinity. 
If this is the case, then the profile of the line should vary 
with the orbital phase of the companion. 

Finally, hydrodynamic modeling indicates that a gas 
giant planet in a circumstellar disk can induce an ec- 
centricity as high as 0.25 in the inner rim that falls off 
as ~ r~^ iKlev fc DirksenI (|2006( ). The modeling also 
predicts that the semi-major axis of the eccentric disk 
precesses slowly (~10°/1000 orbits). Thus the OH line 
profiles should show minimal variation over the orbital 
period of the planet. To determine whether an eccentric 
inner rim can account for the observed asymmetry, we 
model the hot band CO and OH emission line profiles 
assuming that the inner wall has a nonzero eccentricity 
while the outer disk is circular. 



4. MODELING 

In our model, we ad o pt th e disk geometry de- 
scribed bv iBrittain et al.l (I2009D 



puffed up ( DuUemond et a 
of 3.5AU (iBouwrnan et al 
iBenistv et al.l l2010[) 



The 



inner rim is 



.20011 ) with a scale height 
20031 IBrittain et aD [200l 



We adopt the inner radius found 
from the synthesis of the ro-vibrational CO emission 
(13 AU; iBiittain et al. .2003) and the same radial pro- 
file for the disk emission found from our previous mod- 
eling (^R-^/^). We use a stellar mass of 2.4 Mq 
(jvan den Ancker et al.l[T997[ ) , and an inclination of 45°, a 
value within the bounds me asured bv lQuaiiz 
(47° ±2.7°) and lArdila et al. (2Q07) (42° ±5°). 

The magnitude of the bulk velocity of a parcel of gas 
in an eccentric orbit is not axisymmetric and is given by 
v{e)^ = GAh{l + e2 + 2ecos{e))/{a{l - e^)), where G 
is the universal gravitational constant. Mi, is the stellar 
mass, a is the semi-major axis of the eccentric annulus, 
e is the eccentricity, and is the phase of the orbit mea- 
sured from periastron. For non-zero eccentricities, the 
velocity of the gas at periastron is larger than at apas- 
tron, so the side of the line profile arising from gas near 
periastron will extend to higher velocities (Fig. 9). Ad- 
ditionally, the gas in the disk near periastron is closer to 
the star making that gas warmer and brighter. We as- 
sume the temperature scales as T = To{R/ Rq)" and the 
luminosity of the line is proportional to Bi){l — e~'^"). 
Fig. [9] shows the effect of increasing eccentricity on a 
line profile. The synthetic emission line profiles shown 
arise from an eccentric, unresolved annulus inclined by 
45° with a = -1. 

In modeling the line profiles, we assume the line of 
sight falls along the semi-minor axis of the disk and the 
slit is aligned with the semi-major axis of the diskQ Such 
an orientation maximizes the difference between the pro- 
jected velocity on the red-shifted and blue-shifted sides 
of the disk. If the line of sight falls along the semi-major 

^ The position angle of a disk is defined by the apparent semi- 
major axis of the inclined disk. The apparent semi-major axis 
is not generally aligned with the intrinsic semi-major axis of the 
eccentric disk. In this paper the semi-major and semi-minor axes 
refer to intrinsic axes of the disk not the apparent axes. 



axis, this results in the gas at apastron and periastron 
having a zero projected velocity and a symmetric line. 
Thus the eccentricity inferred from our modeling is a 
lower limit. 

For the line profile calculation, the disk is divided 
into annuli, and each annulus is then further divided az- 
imuthally into zones, with each representing a change 
in the projected velocity (along the line of sight) of 
1 km s^^. We assume that the inner wall has an ec- 
centricity that is left as a free parameter. The annuli 
comprising the surface of the disk are assumed to have 
an eccentricity of zero. In the expression for the radial 
temperature variation, T = Tq{R/ Rq)"" , the fiducial tem- 
perature To is assumed to be 1400 K based on what was 
observed for CO. The contribution to the luminosity of 
the emission lines from each zone is proportional to the 
area of the zone and i?£>(l — e~'^"). The emergent spec- 
trum is then calculated by looping over all the zones and 
shifting the luminosity arising from each zone by the pro- 
jected velocity. The spectrum is then convolved with the 
instrument profile, scaled to the flux of the OH feature, 
and resampled into a common velocity space. We leave 
the radial dependence on temperature (and thus the lu- 
minosity) of the OH emission in the inner rim a free 
parameter. 

The fraction of the emission that arises from the in- 
ner wall of the disk is also left as a free parameter for 
the OH and CO distributions. We achieve the best fit 
if 75% of the OH arises from the inner edge, which is 
similar to what has been found for CH+ in HD 100546 
(iThi et al.l[20Tll) . The remaining flux comes from the 
outer disk. For the best fit to the CO, 35% of the lumi- 
nosity arises from the wall and 65 % from the surface o f 
the disk. Fluorescence modeling in IBrittain et al.l (|2009f ) 
found a slightly smaller percentage (15%) of the luminos- 
ity arises from the inner wall. This ratio can be changed 
in their model by increasing the turbulent velocity of the 
gas in the inner rim and decreasing the fiaring or turbu- 
lent velocity of the outer disk. 

We are able to fit the line profiles with the following 
values (and associated 99.7% confidence levels): e — 0.18 
(0.07 to 0.30), a = -2.5 (< -1.2), an OH fraction (aris- 
ing from the inner wall) of 0.75 (> 0.44), and a CO frac- 
tion (arising from the inner wall) of 0.35 (0.19 to 0.55). 
The uncertainty in the eccentricity does not include the 
uncertainty in the viewing geometry. The minimum re- 
duced for fits corresponding to a sampled parameter 
space near these values is 1.0, which corresponds to the 
overall reduced of the model fits shown in Fig. [10] 
& Fig. [TT] If the line of sight does not fall along the 
semi-minor axis of the disk, then the eccentricity of the 
disk will be higher. If the line of sight to the disk is offset 
from the semi- minor axis by 10°, 20°, and 30°, the ec- 
centricity that results in the best fit for each line of sight 
is 0.28 , 0.45, and 0.48, respectively. 

The transition from the eccentric inner wall to the 
circular outer disk is dynamically complicated, but 
models of such disks indicate the transi tion is sharp 
(|Klev fc DirksenI l2006t iRegaly et al.l IMot) . Gas emit- 
ting from the surface of the disk within a few AU of the 
inner rim makes a minimal contribution to the line pro- 
file. As a result, the lack of smooth continuity between 
the eccentric inner annulus and the rest of the circular 
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disk docs not appreciably affect our line profile. 

5. DISCUSSION 
5.1. Massive Companion 

One possibility for the origin of the eccentricity of 
the inner disk is tidal interactions with an embedded 
gas giant planet. When a gas giant planet forms in 
a disk, tidal interactions between the planet and the 
disk caus e the inn er edge of the oute r disk to grow 
eccentric (iKlev fc Di rkscn 2006; Papalo izou et al.ll200ll : 
lLubowlll991[ ). In these hydrodynamic simulations, the 
1:3 Lindblad resonance induces the eccentricity. The 1:2 
Lindblad resonance has a damping effect on the gener- 
ation of the eccentric structure. Thus the gas must be 
cleared from this region of the disk. 

For companion masses greater than q=0.003, where q 
is the ratio of the companion to stellar mass, and for a 
viscosity parameter value of a ^.004, which is typical o f 
circumstellar disks (|Klevl 119990 . iKlev fc Dirksenl (|2006f ) 
show that significant eccentricities (e < 0.25) can be in- 
duced in the portion of the disk exterior to the orbit of 
the embedded gas giant. The azimuthally averaged ec- 
centricity falls off approximately as where r is the 
distance from the star. As a result, only a narrow an- 
nulus of the disk has significant eccentricity (the precise 
radial dependence of the eccentricity depends on the vis- 
cosity of the gas) . 

The possibility of a massive companion to HD 100546 
has been explored previously by several authors. Firstly, 
the transitional nature of the SED may be indicative of 
dyna mical sculpting of the d i sk by a massive con i paniq n 
(e.g. iSkrutskie et al I (Il990l) : iMarsh fc MahonevI lT99l: 



iBrvden et all (I1999D: I Calvet et al.' ('2002^: 'Rice ct all 



(I2003D: iQui-llen et al.l (I2004D: iD'Alcssio c t al. (2005); 
Calvet et al.l (1 2005): ' Furlan et alJ (j2006h . Indeed, 



Kraus fc Ireland (20111) have presented intriguing im- 



agery of the transitional object Lk Ca 15 consistent 
with the presence of a forming planet. iBouwman et al.l 
(|2003l ) model the SED of HD 100546 and conclude 
that the inner ^10 AU of the disk is cleared. These 
authors suggest that the companion must be at least 
5.6 M j in order to open the hole in the disk. Sec- 
ondly, iGradv et al.l (|2005[ ) used long slit images acquired 
with the Space Telescope Imaging Spectrograph on the 
Hubble Space Telescope and confirm that the inner edge 
of the disk is at 13 ±3AU. They also find that the star 
is not centered in the inner hole of the disk. Rather the 
star falls 5 ±3 AU to the northwest. If the offset is due 
to the eccentricity of the disk, then e = 0.38 ±0.24. This 
eccentricity is consistent with our result from modeling 
the line profile of the OH emission. In Paper II, wc dis- 
cuss how the observed variability in the CO fundamental 
emission may be an additional signpost of a forming gas 
giant planet orbiting HD 100546. 

The result of our spectral synthesis of the asymmetric 
OH line profile is s imilar in spirit to the predictions of 
iRegalv et al.l (|2010[ ). although the observations differ in 
detail from their model predictions. Jlcgaly ct al. (2010) 
carried out hydrodynamic simulations of a disk with an 
embedded planet and used the ro-vibrational CO emis- 
sion as a tracer in exploring how disk emission line pro- 
files are modified by the presence of the planet. They 
confirmed the finding of iKlev fc Dirkscn, (20061 ) that an 



embedded planet induces a significant eccentricity in the 
disk. Planets above a critical mass (>3 Mj for a 1 Mq 
star) will cause the outer edge of the gap opened by the 
planet to be elliptical in shape and induce a measurable 
asymmetry in the CO line profiles. Based on our line 
profile modeling, we anticipate detecting a small but ob- 
servable velocity shift in the line profiles due to the dif- 
ference between the projected velocities at apastron and 

periastroii; 

In the IRegalv et all ()2010[ ) predictions for CO emis- 
sion, the disk extends from ~0.2 AU to 5 AU, and CO 
emits over this region of the disk, i.e., within, at, and be- 
yond the orbit of the planet. In comparison, we do not 
measure any emission from the inner optically thin re- 
gion of the disk around HD 100546 (i.e., within 13 AU). 
Instead, the CO lines excited by UV fluorescence arise 
from the eccentric inner rim and the circular outer disk 
extending from ^13-100 AU. The OH emission is mostly 
restricted to the inner rim of the disk. Thus the CO and 
OH line profiles of HD 100546 appear different from the 
synthetic profiles generated bv IRegalv et al.l (|2010[ ). but 
several of the elements of their scenario are consistent 
with the observations. 

5.2. Comparison of OH and [01] Line Profiles 

Our results complement other stud ies of OH disk emis - 
sion in the literature. In their study, iFedele et al.l (|2011[ ) 
detected only the P4.5 line and therefore could not de- 
termine a temperature or the excitation mechanism with 
the available data. The authors do show that the [O I] 
A6300A line and OH emission lines observed toward 
HD 259431 have the same profile, providing some sup- 
port to the hypothesis that the [O I] A6300A line arises 
from the photodestruction of OH. They also find that 
OH emission is more common from flared disks, which 
they suggest indicates the im portance of fluorescen t ex- 
citation. CO observations by iBrittain et al.l ()2007l ) fur- 
ther show the width of the P30 ro-vibrational line to be 
similar to that observed in the OH and [O I] line for 
H D 295431. 

iFedele et al.l ()2011f) also observed asymmetric ro- 
vibrational OH lines toward V380 Ori, which they sug- 
gest could be caused by a deviation from axisymmct- 
ric Keplerian rotation or b y a non-hom ogenous dis- 
tribution of gas. Curiously, lAlecian et al . (2009) find 
that V380 is a spectroscopic binary with a separa- 
tion of rsin(i)<0.33 AU. Asymmetric ro-vibrational CO 
lines similar to the OH ro-vibrational lines presented 
here were observe d in the recent outburst of EX Lupi 
(jCoto et all 120111 ). It should be noted that the EX 
Lupi system is highly variable due to rapid accretion; 
M - 10"'^ Mq yr-^, (|Aspin et al.ll2010( ). Direct com- 
parison to line formation regions in HD 100546 is prob- 
lematic due to HD 100546's comparably slow accretion 
rate; M ~ lO'^ Mq yr'^ , (jCradv et al.ll2005l ). 

The [O I] A6300A emission from HD 100546 has been 
interpreted as arising from the ph otodissociation of OH 
Acke fc van den Ancked |2006[ ). iStorzer fc HoUenbachI 



1998[ ) show that in H/H2 photodissociation regions far 



ultraviolet photons dissociate the OH molecule leaving 
just over half of the oxygen atoms in the ^D2 state - the 
upper level of the [O I] A6300A emission line. The first 
pre-dissociation band of CO is near llOOAwhere the UV 
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field of HD 100546 is quite weak and atomic fiydrogen 
can provide some sliielding. In contrast, tlie first pre- 
dissociation band of OH is centered near I6OOA (l^E-; 
|win Dislr oeck fc Dalgarnolll98^ . wfiicli falls near the 4tli 
positive syste m of CO. Thus CO and OH ca n self-shield 
one another. lAcke fc van den Anckerl ()2006D model the 
line profile of the O I emission from HD 100546 and show 
the gas is distributed from 0.8 - 100 A U. Building on the 
work bv lStorzer fc HollenbachI (|2000( ). they find that the 
surface density of the OH necessary to give rise to the 
O I emission is, 

SoH(fl) = 3.6 X 102i(i?M[/)-2 ^ (3) 

This interpretation requires the inner disk to be 
curiously rich in OH and very depleted in CO 
(N(OH)/N(CO)>100) contr ary to predictio ns for the 
chemistry of inner disks fe.g. lThi fc Bii3[2005l ). 

Adopting this OH surface density and assuming the 
gas is in LTE, we find that the OH would need to be 
as cool as 450 K or have a low density to fall below 
our detection limit. Our upper limit on the rotational 
temperature of the OH is 1400 K which is comparable 
to the temperatur e of the CO a t the inner rim of the 
disk (1400±400 K: lBrittain et all ()2009f )). It is not clear 
why OH in the inner disk would be significantly cooler 
than CO beyond 13 AU. In principle, it is also possible 
that the gas in the inner disk is below the critical den- 
sity necessary to thermalize the vibrational levels such 
that the effective vibrational temperature is much less 
than 1400 K. However, electronic excitation of the OH 
molecule can also populate the excited vibrational lev- 
els. If the OH is not electronically excited because it is 
shielded from the NUV radiation field, then it is unlikely 
that a significant quantity of OH is photodissociated and 
able to give rise to the observed O I emission. Thus it 
seems more likely that the inner disk is largely devoid 
of OH as well as CO and the [O I] A6300A line likely 
arises from a process other than the photodestruction of 
OH. Detailed modeling of the excitation of the OH will 
be considered in future work. 

6. CONCLUSIONS 

We have shown that the observed asymmetric emission 
line profile of OH can be reproduced as arising from an 
inner wall with an eccentricity that is in line with current 
estimates of the eccentricity induced by an orbiting gas 
giant companion. While we cannot conclusively demon- 
strate that the emission arises from an eccentric rim, our 



hypothesis is testable. If the asymmetry is due to the 
eccentricity of the inner rim induced by a massive com- 
panion, then the line profile should remain approximately 
constant in time. One caveat is that if the orbiting planet 
warms the inner rim of the disk in its vicinity, the emis- 
sion line profiles may shift position in velocity space as 
the planet moves through its orbit. If confirmed, this 
observation provides additional circumstantial evidence 
that the gap in the disk at 13 AU has been opened by a 
massive companion. 

Other possible explanations for the line asymmetry in- 
clude an origin in a wind or a turbulent disk atmosphere. 
If the asymmetry is a consequence of disk turbulence, 
the profile shape should vary stochastically with time, 
and we would not continue to see the profile that we now 
observe. If the profile arises in a wind, the asymmetry 
should also be present in other molecular line diagnostics 
that probe similar physical conditions. 

We have also shown that the OH emission lines are sig- 
nificantly narrower than the [O I] A6300A emission line. 
Thus we conclude that the OH and O I arc unlikely to be 
co-spatial. The lack of OH emission from the inner disk 
suggests that the destruction of OH is not responsible for 
the [O I] A6300A emission line, which motivates further 
study into the nature and origin of the emission. 

Based on observations obtained at the Gemini Obser- 
vatory, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc., under a cooperative 
agreement with the NSF on behalf of the Gemini partner- 
ship: the National Science Foundation (United States), 
the Particle Physics and Astronomy Research Coun- 
cil (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research 
Council (Austrafia), CNPq (Brazil) and CONICET (Ar- 
gentina). The Phoenix infrared spectrograph was devel- 
oped and is operated by the National Optical Astron- 
omy Observatory. The Phoenix spectra were obtained 
as part of programs GS-2010B-C-2. We would like to 
thank Bram Acke for providing the published spectrum 
of O I for comparison with our OH spectrum. S.D.B. 
and J.P.L. acknowledge support for this work from the 
National Science Foundation under grant number AST- 
0708899 and NASA Origins of Solar Systems under grant 
number NNX08AH90G. Basic research in infrared as- 
tronomy at the Naval Research Laboratory is supported 
by 6.1 base funding. 



APPENDIX 

APPENDIX A: SPECTROSCOPIC NOTATION 

The total rotational angular momentum is given by the number following the P, the numerical subscript 
corresponds to electron spin of either —1/2 (1) or +1/2 (2), and the +/— signs correspond to even/odd parity. 

Unlike common diatomic molecules such as H2 and CO, the total orbital angular moment of the ground 
electronic state of the OH radical is A=l and is thus doubly-degenerate. The additional degrees of freedom 
require a somewhat more e xtensive labeling scheme to fully describe the transitions. In this paper, we adopt the 
nomenclature described by iSchleicher fc AhearnI ()1982D . Due to the interaction of the electron spin and orbital 
angular momentum (spin-orbit coupling), the ground electronic state is split into "^^(^-1/2 £md '^^^+1/2- The 
total angular momentum, J. is given by the sum of the molecular rotation, = 0,1,2,3,... and electron spin 
s = ±1/2, so that J = -h 1/2 in the X^Hg/j state and J = - 1/2 in the A^Hi/a state. In addition, 
each rotational level is split into two states due to the coupling of the molecular rotation to the orbital motion 
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of the electrons (lambda-type doubling); these split levels are labeled as (+) or (— ) based on the parity of the 
level. The selection rules for a ro-vibrational transition in the ground electronic state are A J = 0,±1 and ± => ^. 

Rotational transitions of J" = J' — 1, J" = J', and J" = J' + 1 are labeled R, Q, and P respectively. For transitions 
with + 1/2 (i.e. the ^113/2 ladder), the rotational designation is labeled with the subscript 2. For — 1/2 the 
rotational designation is labeled with the subscript 1. When the molecule undergoes a transition between the ^113/2 
ladder to the ^lli/2 ladder, the 1 and 2 are given in the order of the transition (upper level first). The parity of the 
ground state is then labeled with the superscript + (for even parity) or — (for odd parity). The ground state electronic 
level is labeled X, the first excited bound electronic level A, the second excited bound electronic level B, etc. In this 
paper wc are only concerned with ro-vibrational transitions within the ground electronic state and thus leave off the 
X. 
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TABLE 1 
Log of Observations 



Date 


Telescope/Instrument 


Target 


Spectral Grasp 


Integration 


S/N 


Airmass 








cm^^ 


minutes 






2010 Dec 22 


Gemini S. /PHOENIX 


Theta Car 


3138 - 3150 


8 




1.5 


2010 Dec 22 


Gemini S. /PHOENIX 


HD 100546 


3138 - 3150 


40 


180 


1.6 


2010 Dec 23 


Gemini S. /PHOENIX 


Theta Car 


2837 - 2850 


8 




1.7 


2010 Dec 23 


Gemini S. /PHOENIX 


Theta Car 


2027 - 2037.5 


8 




1.4 


2010 Dec 23 


Gemini S. /PHOENIX 


HD 100546 


2837 - 2850 


28 


180 


1.8 


2010 Dec 23 


Gemini S. /PHOENIX 


HD 100546 


2027 - 2037.5 


32 


79 


1.4 



TABLE 2 
Equivalent Widths of OH lines 



Line Wavenumber Equivalent Width 



cm ^ 10 •^cm ^ 



P9.52 - 3146.18 3.8±0.5 

P9.52+ 3145.49 4.1±0.5 

P10.5i+ 3142.06 3.8±0.5 

P10.5i- 3141.04 5.3±0.5 

P15.52- 2840.82 <.5 

P15.52+ 2839.58 <.5 

P16.5i+ 2839.03 <.5 

P16.5i- 2837.43 <.5 



TABLE 3 
Unblended Lines 



Line 


Wavenumber 




cm^^ 


v=3-2 P15 


2030.16 


v=6-5 R5 


2032.83 


v=4-3 P8 


2033.14 


v=l-0 P16 "CO 


2033.42 


v=3-2 P14 


2034.41 
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Fig. 1. — High resolution near infrared spectra of HD 100546 in the observed frame. The top panel shows the ratioed M— band spectra 
containing the ro-vibrational CO transitions. Numerous emission features are evident and are labeled. The positions of these features in 
this spectrum are marked by vertical dot-dashed lines. The bottom panel shows the observed spectrum for this wavelength range and the 
telluric standard in red. The strong emission visible in the top panel is also clearly seen here. The gaps in the spectrum are of regions 
where the transmittance of the atmosphere is less than 50%. The fundamental ro-vibrational CO P26 line and numerous hot band features 
are evident. 
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Fig. 2. — High resolution near infrared spectra of HD 100546 in the observed frame. The top panel shows the L— band spectra containing 
the ro- vibrational OH transitions while the bottom panel presents the observed spectrum (black) and the telluric standard (red). The gaps 
in the spectrum are of regions where the transmittance of the atmosphere is less than 50%. The positions of the v=l-0 OH lines are labeled 
with vertical dot-dashed lines. The P15.52± and P16.5i± ro-vibrational emission lines are not detected. We use the noise in the continuum 
to set an upper limit on the equivalent widths of these lines, which are then used to put an upper limit on the rotational temperature. 
Detections of emission lines in this range arc below the 1-cr level. 
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Fig. 3. — High resolution near infrared spectra of HD 100546 in the observed frame. The top panel shows the ratioed L— band spectrum 
containing the ro- vibrational OH transitions while the bottom panel presents the observed spectrum (black) and the telluric standard (red). 
The gaps in the spectrum are of regions where the transmittance of the atmosphere is less than 50%. The positions of the v=l-D OH lines 
are labeled with vertical dot-dashed lines. The P9.52± and P10.5i± ro-vibrational emission lines are detected and labeled with dot-dashed 
vertical lines. 
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Fig. 4. — Excitation diagram of ro-vibrational OH emission lines. The quantity ln(F/hc£'Ag) is plotted against E/k, the energy of the 
upper level, where F is the line flux (in erg/s/cm^/cm~^), i> is the wavenumber of the transition (in cm~^), g is the multiplicity of the 
upper level, and A is the Einstein A cocfhcient of the transition. Thus the rotational temperature of the gas is given by the negative 
reciprocal of the slope. The solid lino shows the fit to the detected values and the upper limits. The Icr upper limit on the rotational 
temperature of the OH gas, based on the best fit line, is 1400 K. 
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Fig. 5. — Individual OH lines used to construct the average OH line profile in the rest frame of the star. The S/N of the average 
continuum is ~325. The average OH line to continuum ratio is ~ 3% and the continuum noise is ~ 0.6%. The HWZI is 12 ±2 km s~^. 
The asymmetric line shape is evident in most of the individual lines, with approximately 80% of the emission in the blue and 20% in the 
red. The P10.5j^+ line was shifted 1.5 km s^'^ blueward to line up with the other OH lines such that residuals between the P10.5j+ line 
and the other OH lines corresponded to the noise along the noise along the continuum. 
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Fig. 6. — Unblended hot band CO lines in the rest frame of the star. The signal to noise of the continuum of the average profile is ~110. 
The HWZI of the line is 13 ±2 km . Line positions are reported in Table 3. 




Velocity (km s ^) 

Fig. 7. — Comparison of the [O I] A6300A line and the average OH line profile (adapted from lXcke et all l|2005l ')'l in the rest frame of 
the star. The scale of the [O I] A6300A line (black) is the left vertical axis. The scale of the OH emission line (red) is denoted on the right 
vertical axis. The line profiles are clearly different, and the HWZI velocity of the OH emission is much narrower than that of the O I line 
indicating that it does not arise from the same region of the disk as the OH emission. 




Fig. 8. — Comparison of the average hot band CO hne profile (black) and the average OH line profile (red) in the rest frame of the 
star. The scale of each emission line profile is plotted in the loft axis. The line profiles are clearly different, but the width of each line is 
consistent. 
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Fig. 9. — Synthetic, spatially unresolved line profiles of emission from an eccentric inner wall that is inclined by 45°to the line-of-sight in 
the rest frame of the star. The other model parameters are alpha =-1, an radius for the inner wall of 13 AU, and a stellar mass of 2.4Mq 
(see text for details). The emission is entirely from the wall component in this example to best illustrate the asymmetry, which grows as 
we increase the eccentricity of the wall. For this example the equivalent width is held constant and equal to the measured equivalent width 
of the average OH line. The asymmetric emission bump is blueward of line-center due to the direction we have chosen for the rotation of 
the gas and the viewing angle of disk; these choices are entirely general. For HD 100546 we see a similar blueward emission bump in the 
OH omission (see text). 
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Fig. 10. — Synthesized spectrum of OH emission arising from the inner wall of a disk with an eccentricity of 0.18 in the rest frame of the 
star (red line). The average OH line profile (black line) is plotted in the upper panel. In the lower panel we plot the difference between the 
data and the model. To reproduce these line shapes we assumed that the ratio of the emission luminosity of the wall to the luminosity of 
the disk was approximately 3:1. Lower ratios produce a more double-peaked structure, while higher ratios do not adequately fill out the 
red component of the OH emission line. The minimum reduced chi-squarcd value, /v, is 1.0. 
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Fig. 11. — Synthesized spectrum of CO emission arising from the inner wall and the disk in the rest frame of the star (red line). The 
average line profile of the unblended hotband CO emission lines (black line) is plotted in the upper panel. The difference between the 
observed line profile and the modeled line profile is plotted in the lower panel. The synthetic CO line profile includes components from both 
the inner wall and disk. We assumed that the ratio of the emission luminosity of the wall to the luminosity of the disk was approximately 
1:3. The reduced chi-squarcd value, X'^/^: is l-l- 



